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(57) ABSTRACT 

Method and system for analyzing and processing information 
on one or more aircraft flight paths, using a four-dimensional 
coordinate system including three Cartesian or equivalent 
coordinates (x, y, z) and a fourth coordinate 6 that corre- 
sponds to a distance estimated along a reference flight path to 
a nearest reference path location corresponding to a present 
location of the aircraft. Use of the coordinate 8, rather than 
elapsed time t, avoids coupling of along-track error into air- 
craft altitude and reduces effects of errors on an aircraft 
landing site. Along-track, cross-track and/or altitude errors 
are estimated and compared with a permitted error bounding 
space surrounding the reference flight path. 
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Climb Altitude Profile 
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aircraft type: B727 
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Along-Track Position during Descent 
aircraft type: B727 
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TRAJECTORY SPECIFICATION FOR HIGH 
CAPACITY AIR TRAFFIC CONTROL 

ORIGIN OF THE INVENTION 

5 

This invention was made, in part, by one or more employ- 
ees of the U.S. government. The U.S. government has the 
right to make, use and/or sell the invention described herein 
without payment of compensation therefor, including but not 
limited to payment of royalties. 10 

CROSS REFERENCE TO RELATED 
APPLICATIONS 

FIELD OF THE INVENTION 15 

This invention relates to control of aircraft traffic through 
specification of trajectory coordinates. 

BACKGROUND OF THE INVENTION 20 

Use of a four-dimensional description of an aircraft trajec- 
tory, involving three Cartesian coordinates (x, y, z) and 
elapsed time t, was proposed in 1972 by H. Erzberger and T. 
Pecsvaradi (“4D Guidance System Design With Application 
To STOL Air Traffic Control”, Yb th Joint Automatic Control 
Confer., Stanford, Calif. Aug. 16-18, 1972, pp. 445-442) and 
has been considered by many other workers since then. One 
potential problem with use of (flight) time as a fourth coor- 3Q 
dinate is that an along-track error couples into aircraft alti- 
tude. In this approach, when an elapsed time specified by 
flight plan coordinates is reached or exceeded, an aircraft may 
be required to land many miles short of, or many miles 
beyond, its destination. Alternatively, failure of an air traffic 35 
management (ATM) ground computer could abruptly deposit 
responsibility for safe separation of aircraft proceeding along 
the same general route onto human air traffic controllers 
(ATCs). One ultimate goal of an automated ATM system is to 
remove, or minimize reliance on, ATCs from active partici- 4Q 
pation in maintenance of separation between aircraft. 

What is needed is a method and system for specification of 
an aircraft trajectory in terms of coordinates that (1 ) allow use 
of higher capacity of traffic in a given volume of airspace, (2) 
define bounds on allowable along-track, cross-track and ver- 45 
tical errors, and (3) compensate for deviations from a time 
schedule so that an aircraft is not required to fly with unreal- 
istic velocity or unrealistic angle parameters and is not 
required to execute a landing procedure at a location that is 
spaced far apart from a destination location. The approach 50 
should allow for in-flight changes in flight parameters to take 
account of a changed environment. 

SUMMARY OF THE INVENTION 

55 

These needs are met by the invention, which provides a 
four-dimensional location coordinate system, for example, 

(x, y, z, 6), where 6 is not elapsed time but is a measured or 
calculated distance along an aircraft flight path. Along-track 
error, cross-track error and altitude error are calculated rela- 60 
five to an ideal, conflict-free reference flight path. Elapsed 
time is monitored and used as an independent variable for 
computing or measuring along-track, cross-track and altitude 
locations and errors. An error tube is numerically constructed 
around a central axis, such as the conflict-free reference flight 65 
path, with an along-track range AAT, a cross-track range ACT 
and an altitude range AAL that may vary with the distance 5 


2 

along the reference flight path. Altitude is preferably specified 
as a function of along-track distance 8, rather than as a func- 
tion of time t. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a curvilinear segment of, and grid for, a 
flight path during an aircraft turn. 

FIG. 2 illustrates a method for determining a location on a 
reference flight path corresponding to an aircraft present loca- 
tion. 

FIG. 3 illustrates determination of reference altitude and of 
actual altitude as a function of along-track (AT) location. 

FIGS. 4 and 5 graphically illustrate synthesized altitude 
profile and a best-fit parabola and error bounds for a constant 
calculated air speed (CAS) for an ascending aircraft (Boeing 
757) 

FIG. 6 graphically illustrates synthesized altitude profile 
and a best-fit parabola and error bounds for a constant calcu- 
lated air speed (CAS) for descending aircraft (Boeing 727). 

FIG. 7 graphically illustrates synthesized altitude profile 
and a best-fit parabola and error bounds for a constant calcu- 
lated air speed (CAS) for an ascending aircraft (Boeing 727) 

FIG. 8 graphically illustrates an aircraft turn geometry. 

FIG. 9 graphically illustrates error bounds for an aircraft 
that is transitioning from ascending to level flight. 

FIG. 10 graphically illustrates an aircraft air speed enve- 
lope versus pressure altitude for an MD-80 aircraft. 

FIGS. 11, 12 and 13 graphically illustrate cockpit displays 
of traffic information (CDTIs) for horizontal and vertical 
guidance and for situational awareness for a representative 
aircraft flight. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

A trajectory specification used in the approach disclosed 
here will be used for communicating trajectory information 
between an aircraft and one or more ground systems, such as 
an Airline Operation Center (AOC), for scheduling, maintain- 
ing aircraft separation and other relevant activities. The tra- 
jectory specification has several requirements. The specifica- 
tion: 

(i) should be able to precisely specify any “reasonable” 4D 
assigned trajectory, 

(ii) should be able to precisely specify error tolerances 
relative to the assigned trajectory, 

(iii) should be based on a global earth-fixed coordinate 
system, 

(iv) should be parametric and reasonably compact, 

(v) should be based on a text format readable by humans, 
and 

(vi) should be suitable for an international standard. 

The first requirement is that the format be able to precisely 
specify any “reasonable” 4D trajectory (3D location as a 
function of time). A unique 3D location must be precisely 
determined at each point in time, and the set of specifiable 
trajectories must not be unreasonably restrictive. Efficient 
climbs and descents must be allowed, for example, and turns 
must be allowed during climb and descent. The horizontal 
path will be restricted to straight (great circle) segments con- 
nected by turns of constant radius to simplify computations 
and conformance monitoring. These restrictions should not 
significantly limit practical routing flexibility. Note that wind 
optimal routes can be approximated with sufficient accuracy 
for practical purposes using great circle segments of, say, 100 
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to 200 nautical miles (mni) in length, depending on the length 
of the flight. More general horizontal path segment types can 
be added later if desired. 

The format should also have an ability to specify error 
tolerances for the flight technical error in each of the three 
axes: along-track (AT), cross-track (CT), and altitude or ver- 
tical (ALT). The error tolerances will precisely determine a 
3D bounding space in which the aircraft is required to be 
contained at any point in time. Those bounds will be the key 
to assuring that the minimum required separation is main- 
tained at all times without the attention of a human controller. 
If an aircraft fails to conform, or is expected to fail shortly its 
status will be temporarily downgraded to unequipped; and if 
necessary it will be automatically issued a basic heading or 
altitude resolution advisory (by a system that is outside the 
scope of this disclosure). After the situation is under control, 
the aircraft will either reacquire a new assigned trajectory or 
be handed off to a human controller for conventional separa- 
tion monitoring until it is able to reacquire a new assigned 
trajectory. 

Trajectories will be synthesized to guarantee the minimum 
required separation for a specified period of time called the 
conflict time horizon, which could be perhaps fifteen minutes. 
A key point is that, if the trajectories are correctly synthe- 
sized, conformance by any two aircraft will guarantee the 
minimum required separation distance between the aircraft 
for a specified period of time, regardless of where each air- 
craft is within its own bounding space. Otherwise stated, the 
bounding spaces themselves will always maintain the mini- 
mum required separation. Note that minimum separation 
standards are specified ir terms of the separation distance 
between aircraft, regardless of velocities or higher-order 
dynamics. Hence, the trajectory error tolerances will also be 
specified in terms of distance or length, as a measurement 
unit. Velocity and acceleration can affect future conformance, 
but actual present conformance will not depend on these 
quantities. Nevertheless, a conformance monitoring system is 
free to use velocity and acceleration to attempt to predict 
imminent nonconformance. 

In the current air traffic system, standard navigational con- 
formance bounds of ±4 nautical miles (mni) in cross-track 
(CT) defines a lane width of 8 mni. However, these bounds are 
routinely violated for various reasons, such as loose piloting 
or controllers issuing heading “vectors” but not entering them 
into the system. In the vertical or altitude (ALT) axis, con- 
formance bounds apply only in level flight, and no bounds 
apply in the along-track (AT) axis, except arrival time con- 
straints. The lack of rigorous conformance bounds in a con- 
ventional system makes conformance monitoring a “fuzzy” 
problem, which Reynolds and Hansman (“Investigating Con- 
formance Monitoring Issues In Air Traffic Control”, ICAT- 
2003-5, Inti. Center for Air Transportation, M.I.T., November 
2003) have attempted to solve using fault detection methods. 
However, conformance monitoring is precisely defined if 
conformance bounds are based on location only and are 
specified precisely. The more difficult and “fuzzy” problem is 
the detection of faults that could lead to imminent non-con- 
formance, and that is where Reynolds and Hansman approach 
may still apply. 

The error tolerances will be based on Required Navigation 
Performance (RNP) specifications, but these tolerances could 
be relaxed in sparse traffic . Because wind velocities cannot be 
modeled exactly, the most challenging axis for which to set 
tolerances is the AT axis. Tightening the AT tolerance 
increases airspace capacity, but it also increases the probabil- 
ity that aircraft will be required to fly at inefficient or unreal- 
istic airspeeds. Along-track location error tolerances must be 
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set as a compromise between those two effects. For more 
flexibility, AT location error will be allowed to vary linearly 
with time. Also, the assigned (reference) AT location and 
velocity will be updated periodically to compensate for errors 
5 in modeling and prediction of AT wind magnitudes, but only 
when doing so does not result in a conflict. The management 
of AT location assignments and tolerances wifi be discussed 
in more detail in the following. 

A third requirement is that the format be based on a global, 
1 0 earth-fixed coordinate system, which will provide a common 
reference. Local coordinate systems, such as the (pseudo- 
Cartesian) stereographic projection used within each Air 
Route Traffic Control Center (ARTCC), are inappropriate for 
enroute monitoring because these coordinate systems are 
1 5 each valid only within one Center. The complexity of switch- 
ing coordinate systems for each Center would be unnecessar- 
ily complicated. The standard WGS84 geodetic coordinate 
system (latitude, longitude, and altitude) will be used as the 
reference coordinate system for enroute airspace. Use of a 
2 0 local coordinate system may be convenient in a terminal area, 
however, and this option will also be available. Also, a curvi- 
linear flight path coordinate system is introduced for speci- 
fying and monitoring the flight technical error tolerances. 

25 A fourth item in the requirements list is that the format be 
parametric and reasonably compact. A continuous 4D trajec- 
tory can be approximated by a simple sequence of discrete 4D 
points (x, y, z, t), but such an approach tends to be inefficient 
in terms of storage and use of bandwidth. This discrete 
30 approach also fails to capture the structure of the trajectory. 
Real trajectories consist of discrete segment types, such as 
climb at constant CAS (Calibrated Airspeed), cruise at con- 
stant Mach number, etc., but discrete 4D points do not convey 
that structure. This disclosure uses a structured, parametric 
35 approach based on straight (great circle) segments, constant- 
radius turn segments, and low-order polynomial approxima- 
tion. In addition to data compression, polynomials also pro- 
vide other desirable features, such as analytic differentiation 
and simplified interpolation. 

40 A more fundamental problem with using a sequence of 
discrete 4D points is that AT location error couples into alti- 
tude. Assume, for example, that an aircraft is on approach for 
landing and is one minute ahead of schedule, but still within 
tolerance. If altitude is specified as a function of time, the 
45 aircraft will be required to land several miles before it reaches 
the runway! On the other hand, if altitude is a function of AT 
location, the aircraft will be required to land at the runway 
regardless of its status. Clearly the latter approach is prefer- 
able. Discrete 4D points are suitable for specifying aircraft 
50 trajectories that have already been flown, but this approach is 
the best choice for specifying trajectories, with error toler- 
ances, yet to be flown. 

A fifth requirement is that the format be in plain text, 
readable by humans. The traditional standard for text is 
55 ASCII (American Standard Code for Information Inter- 
change), but the new more general standard is Unicode. For 
purposes of this disclosure, the ASCII subset of Unicode is 
sufficient. Text-based formats typically provide less efficient 
storage than binary formats, but a text -based format is often 
60 more flexible and less subject to error. Also, a text-based 
format is more convenient because it can be read directly by 
humans. Text can be compressed and encrypted into a binary 
format for efficient and secure radio transmission, but must be 
decompressed and decrypted at the receiving end to recover 
65 the original text. XML, the Extensible Markup Language, is 
a new standard text-based format for specifying structured 
data and transferring it across platforms. An XML data format 
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standard can be used for specifying trajectories for appropri- 
ately equipped aircraft in future high-capacity airspace. 

XML is designed for creating application- specific or 
domain- specific standards for data specification and transfer. 
The resulting text-based standards are intended to be inde- 
pendent of any particular computer platform or language. 
XML is rapidly replacing binary formats for automated busi- 
ness -to -business transactions and is being widely used for 
computing standards such as Scalable Vector Graphics 
(SVG). Whereas binary formats typically require the same 
data to be transferred in the same precise order every time, 
XML provides more flexibility in the selection and ordering 
of the data fields. The flexibility of XML will be indispens- 
able for trajectory specification because each trajectory can 
have a variable number of segments of various types. The 
flexibility will also allow trajectories to be updated without 
repeating all the data that remains unchanged from the pre- 
vious update, which could more than compensate for the 
inherent inefficiency of text-based data. Note also that XML 
text compresses well for efficient use of bandwidth. 

The sixth requirement for the proposed trajectory specifi- 
cation standard is that it be suitable for an international stan- 
dard that is recognized by, and can be automatically flown by, 
any standard flight management system (FMS). The standard 
will be used onboard aircraft to downlink requested trajecto- 
ries constructed by the FMS or constructed by the pilot using 
a graphical user interface. The standard will also be used on 
the ground to check for conflicts and to uplink assigned tra- 
jectories. Developing a consensus for an international stan- 
dard is obviously a major challenge, but such a common 
language can greatly simplify the logistics of high-capacity 
ATM. With a common trajectory language, the chances of 
mi scorn muni cation will be much less than they would be 
without one. If adopted, the actual communication mecha- 
nism would probably be an extension of CPDLC or a new data 
link message over the Aeronautical Telecommunication Net- 
work (ATN): Estimation of Flight Path Errors. 

FIG. 1 illustrates a segment of an aircraft reference (error- 
free) flight path that is followed during a turn, at substantially 
constant altitude, by a turn angle (|), showing AT error bounds, 
CT error bounds, and use of two integrated coordinates sys- 
tems (Cartesian/polar/Cartesian) for a preceding linear seg- 
ment LSI, a turn segment TS, and a succeeding linear seg- 
ment LS2 for this portion of the reference flight path. 

A nearest -point reference present location NRPLA of the 
aircraft can be determined from the actual present location 
APLA, as illustrated in FIG. 2, which also illustrates the 
reference flight path RFP and the actual flight path AFP. A 
perpendicular segment PerpS, relative to the reference flight 
path RFP, is drawn or otherwise constructed from the APLA. 
An intersection of PerpS and RFP, referred to herein as a 
“perpendicular foot,” of the APLA on the reference flight path 
RFP determines the NRPLA. The NRPLA is unique if the 
aircraft is in conformance or if certain conditions are met. 

If the NRPLA is not unique so that M such corresponding 
locations exist (Mi=2), one may use (1) the “earliest” location 
NRPLA that corresponds to the APLA, (2) the “latest” loca- 
tion NRPLA, or (3) a selected linear combination of all cor- 
responding locations RPLA(m) (m=l, . . . , M) that corre- 
spond to the APLA, as set forth inEqs. (1-1), (1-2) and (1-3), 
respectively: 

(x(NRPLA) ANRPLA) ^(NRPLA)=(x(RPLA(w =l),y 

(RPLA(m=ll)/(RPLA(m=l)), (1-1) 


(x(NRPLA) ANRPLA) ^(NRPLA)=(jc(NRPLA(k =N), 

XNRPL A(k =N) ,z(NRPL A(« =N) (1-2) 


6 

(x(NRPLA),^(NRPLA)^(NRPLA)M=2w m (x(NRPLA 

(m),XNRPLA(m),z(NRPLA(m)), (1-3) 

m=l 

5 where w m is a non-negative weight value for the mth NRPLA 
location on the reference flight path RFP. Other prescriptions 
for determining the coordinates (x(NRPLA), y (NRPLA), 
z(NRPLA) corresponding to the actual present location 
APLA. The distance 8 of the NRPLA location along the flight 
1 0 path can be defined as an integral of differential arc length ds 
along the reference flight path RFP from an initial location to 
the NRPLA location: 

6=6(NRPLA)=W (W ^ ) ^. (2) 

15 Another quantity of interest here is a distance, denoted 
5(RPLAET) and measured along the RFP, from the initial 
location to an elapsed time aircraft location, where the aircraft 
would be at the present elapsed time t, if the aircraft had flown 
with zero error along the RFP: 

20 

6 (RP 'LAET)=j RFP ( RPLAET) ds . (3) 

The value, 8=8(NRPLA), is processed and used as a fourth 
coordinate in a description of the actual flight path. 

With reference to FIG. 1, a CT error value 8 cr can be 
estimated by a distance d between the APLA and NRPLA 
locations: 

6 cr =rf{APLA, NRPLA} (4) 

An AT error value b AT can be estimated by an arc length 
3 0 difference between the locations NRPLA and RPLAET, mea- 
sured along the RFP 

b A f= I6(RPLAET)-6(NRPLA) I . (5) 

35 A z-difference, z(RPLAET)-z(APLA), can be interpreted as 
the altitude error: 

b ALT = lz (RPLAET) -z( APLA) I (6) 

The three error components, b ATi b CT , and 8^ r , are consistent 
40 with each other, and involve independent components. 
Because of the definitions of b AT 2 and b CT , the Euclidean 
distance d(APLA, RPLAET) between the actual aircraft loca- 
tion APLA and the elapsed time (error-free) location 
RPLAET will be close to, but may not have the same value as, 
45 the grid distance 

{&at 2 +&ct 2 +&alt} 1/2 - 

The reference flight path RFP is a path the aircraft would 
follow in an ideal situation of zero technical error. Error 
50 tolerances, specified for each of the AT, CT and ALT axes, 
define a tube surrounding the RFP. So-called “straight” seg- 
ments, associated with minimum distance between geodetic 
points are great circles but may be close to linear in latitude 
and longitude for short segments. In a straight (great circle) 
55 segment of a reference trajectory, the local flight path coor- 
dinate system is approximately Cartesian within the range of 
practical error tolerances, and the AT and CT coordinates of a 
point can be determined with established great circle algo- 
rithms. 

60 The curve representing the reference flight path RFP is a 
continuous curve including primarily (1) one or more great 
circle (GC) segments and (2) one or more turn segments with 
substantially constant radius, both referenced to substantially 
level flight and (3) one or more substantially ascent or descent 
65 segments for change of aircraft flight level (FL). The arc 
length integral value for b can be expressed as a sum of two 
components, the first being the accumulated arc length for the 
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consecutive segments already completed. This component is 
straightforward to compute, because the consecutive compo- 
nents are continuous and are joined end-to-end, and each of 
these completed segments is fully characterized by the 
parameters describing the corresponding segment. 

The second component is a not-yet-completed segment, 
the segment that is presently being followed. This segment is 
one of the types ( 1 ) or (2) or (3) characterized in the preceding 
paragraph, and the presently increasing arc length of this 
segment is determinable using the parameters that character- 
ize this segment type. For a GC segment, the differential arc 
length increases according to 

ds { =R A A<$>, (7) 

where A(|), refers to the angle differential subtended at the 
Earth’s center between the beginning of this (incomplete) 
segment and the present location of the aircraft (traveling 
along this GC segment) and R^ is the constant radius of this 
GC segment. For a turn segment, the differential arc length 
increases according to 

ds 2 =R 7 Ar|, (8) 

where Ar| refers to the angle through which the aircraft has 
turned and R r is the constant radius of the turn. For an FL 
change segment, the differential arc length increases accord- 
ing to 

ds 3 =\Akcscy\, (9) 

where y is the substantially constant angle, relative to a local 
horizontal plane, at which the aircraft is ascending or 
descending and Ay is the cumulative altitude increase or 
decrease, relative to the initial point of FL change for this 
segment. Coordinate Systems and Transformations. 

An assigned trajectory consists of a 4D reference trajectory 
and flight technical error tolerances. The reference trajectory 
is the precise 4D trajectory the aircraft would fly in the ideal 
case of zero flight technical error. It is a collection of precise 
3D locations that vary as a function of time, and the location 
at any point in time will be referred to as the reference loca- 
tion. The error tolerances, on the other hand, are the maxi- 
mum allowed error in each of the three axes: AT, CT, and ALT. 
These tolerances define a 3D bounding space around a time- 
based sequence of reference locations, within which the air- 
craft must remain within to be in conformance. 

A WGS84 geodetic coordinate system is preferably used as 
a global standard for specifying reference trajectories. 
Straight (i.e., m inimum distance) segments between geodetic 
points are great circles in general, but for short segments 
(spaced apart from the earth’s poles) a great circle is substan- 
tially linear in latitude and longitude. Geodetic coordinates 
are inconvenient for specifying and monitoring error toler- 
ances, however. For that purpose, a curvilinear flight path 
coordinate system, which follows the reference trajectory, 
will be used. An example of a segment of such a curvilinear 
flight path coordinate system is illustrated in FIG. 1, showing 
the AT and CT grid. 

A curvilinear flight path coordinate system is a combina- 
tion of Cartesian and polar coordinate systems. The first step 
in converting from WGS84 coordinates to curvilinear coor- 
dinates is to determine the type of the local coordinate region, 
which is Cartesian in the (reference) straight segments and 
polar (or cylindrical in 3D) in the assigned turn segments as 
shown in the Figure. Actually, these regions are not strictly 
Cartesian or polar, because the regions follow the curvature of 
the earth; but for practical purposes they are Cartesian or 
polar within the local region of reasonable flight technical 
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errors. The key point is that each segment defines its own local 
coordinate system, which is Cartesian for straight segments 
and polar for turn segments. 

Note that the bounding space, BS in FIG. 1, is based on the 
5 definition of the AT and CT and ALT error coordinates . Thus , 
the bounding space follows the curvature of the flight path and 
may have curved, not planar, surfaces, as indicated in FIG. 1. 
An alternative way to define the error coordinates in terms of 
a time varying Cartesian coordinate system that follows the 
1 0 aircraft as it turns , but that would not work very well . With thi s 
alternative definition, the bounding space shown in FIG. 1, for 
example, would be a rectangular parallelepiped and would 
not conform to the curvature of the flight path. As the refer- 
ence location progresses around the turn, the corresponding 
15 bounding space would be a rectangle that rotates with the 
reference track angle, which would clearly be inappropriate. 

The AT and CT coordinates of two points are defined 
relative to a great circle or relative to a circular turn segment, 
2Q connecting the points. Great circle equations apply only in the 
Cartesian-coordinate regions of the curvilinear flight path 
coordinate system, which correspond to the straight (great 
circle) segments of the assigned trajectory. However, these 
equations can easily be adapted for use in the polar coordinate 
25 regions also, which correspond to the turning segments of the 
assigned trajectory. One approach begins by computing the 
AT and CT coordinates as if a point were still in the preceding 
Cartesian region, then converts to polar coordinates. The ori- 
gin of the polar coordinate system will be the center of the 
30 turn arc, and the reference azimuth angle is determined at the 
start of the turn. The actual CT coordinate will be the radial 
coordinate minus the nominal radius of the turn, so that the 
reference CT coordinate is always zero, consistent with the 
straight segments. The AT coordinate will be the angle rela- 
35 tive to the start of the turn, multiplied by the nominal radius of 
the turn. Note that if the aircraft is flying the turn with a CT 
error, the actual radius of the turn will be different than the 
nominal radius, and the actual AT distance traveled by the 
aircraft will be different than the AT coordinate. 

40 A 4D trajectory also includes a vertical profile, which is 
altitude as a function of time or AT location. While either time 
or AT location could be used as the independent variable AT 
location provides a critical advantage: it fixes the reference 
trajectory in the earth-fixed coordinate system, which simpli- 
45 fies conflict calculations. Using time as the independent vari- 
able, on the other hand, would allow the reference trajectory 
to drift (relative to the earth-fixed coordinate system) with the 
AT location error. As mentioned earlier, a trajectory can be 
visualized as a 3D tube in space, through which the aircraft 
50 flies, with the AT location in the tube being the fourth coor- 
dinate. The tube itself should be fixed in space. Making alti- 
tude a function of time causes the tube to shift in space as a 
function of the AT error. 

The AT location will be specified as a low order polynomial 
55 function of time for each segment, as will be discussed in the 
next section. Also, altitude will be specified as a low-order 
polynomial function of the actual (as opposed to reference) 
AT location, as illustrated in FIG. 3, which shows the refer- 
ence trajectory as a solid curve with a dot on a curve. The 
60 other dot, in the upper right portion of the F igure, indicates the 
actual location of the aircraft. The AT location error is shown 
as a difference between the actual and reference locations. 
Similarly, the ALT error is shown as a difference between the 
actual and reference altitudes. The key point of FIG. 2 is that 
65 the reference ALT is a function of the actual, rather than 
reference, AT location. As noted in the preceding, this fixes 
the 3D flight tube in the earth-fixed coordinate system. 
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In case this distinction is still unclear or seems unneces- 
sary, a simple thought experiment should help clarify and 
justify it. Consider a trajectory that is specified all the way 
through final approach to actual landing. Now suppose the 
flight is one minute behind schedule (but still within toler- 
ance). If altitude is specified only as a function of time, then 
the aircraft will be required to touch down one minute — or 
several miles — before it reaches the runway! That obviously 
won’t work. But if ALT is specified as a function of the scalar 
AT location, the flight will be required to touch down at the 
correct point on the runway regardless of how far ahead of or 
behind schedule the aircraft may be. Clearly, the latter 
approach is correct. 

Polynomial Approximation 

In the current air traffic system, ALT profiles are difficult to 
predict accurately based on information available on the 
ground. Part of the problem is that weight and thrust (or 
throttle setting) are not accurately known on the ground. 
Another major source of altitude prediction uncertainty is the 
lack of knowledge of the actual time of initiation of an altitude 
transition. When an ascend or descend command is received, 
the time taken by the pilot to initiate the maneuver can vary by 
up to about a minute. As a result, controllers must reserve a 
large block of airspace around any aircraft that is in, or is 
about to enter, an altitude transition. With better information 
available on the ground, altitude can be assigned more pre- 
cisely, which will increase airspace capacity. 

An ALT profile is specified to provide reasonable bounds 
on altitude without significantly compromising efficiency. In 
the absence of conflicts, the ALT profile approximates a pro- 
file that the aircraft would be most likely to fly, if uncon- 
strained (or with) normal arrival time constraints. In ascent 
and descent, commercial transport airplanes normally fly 
with the throttle fixed and with feedback to the elevator to 
maintain constant CAS (at lower altitudes) or constant Mach 
number (at higher altitudes). In the future, the intended CAS/ 
Mach schedule will be known on the ground, as will the 
throttle setting and the estimated weight of the aircraft. The 
predicted wind, temperature, and pressure fields will also be 
available from a centralized weather service. Given this data, 
the ALT profile can usually be predicted fairly accurately, and 
an approximation of the predicted ALT profile can be used as 
the assigned ALT profile. If the wind data are reasonably 
accurate, and if the ALT tolerances are reasonable, the aircraft 
should be able to conform to the specified trajectory by flying 
as usual with the specified power and CAS/Mach schedule. 

The mechanics of the procedure will be similar to what is 
currently done in the Center/TRACON Automation System 
(CTAS), but with a few key differences. CTAS is a suite of 
ATC/ATM decision support tools that is being developed at 
NASA Ames Researcl Center. CTAS currently has to guess at 
the weight and the CAS/Mach schedule to be flown, but these 
data will be available from the aircraft or from the AOC. A 
more fundamental difference is that the predicted trajectory 
will actually become the assigned trajectory if the trajectory 
is free of conflicts; otherwise, the trajectory will be modified 
to eliminate any conflicts before becoming the assigned tra- 
jectory. The current ATC system has no such precisely 
defined ALT profiles, and use of a CTAS prediction does not 
guarantee actual, or even attempted, conformance. In fact, the 
notion of vertical conformance itself isn’t even defined for an 
altitude transition. 

The CTAS software process that predicts trajectories is 
referred to as a Trajectory Synthesizer (TS). The TS contains 
performance models of all major aircraft types, and types that 
are not modeled directly are approximated with similar avail- 
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able models. The inputs to the TS for each aircraft include the 
aircraft type and weight, CAS/Mach values, throttle settings, 
the flight plan, and the current weather data file from the 
Rapid Update Cycle (RUC). The output is the predicted 4D 
5 trajectory in the form of a discrete series of coordinates in 
which the time increment varies with the dynamic state. The 
TS or its functional equivalent can be used to construct the 
reference trajectory as closely as possible to the trajectory 
that would have been flown without the constraints, 
to FIG. 4 shows the ALT profile synthesized by the TS for a 
constant-CAS climb segment of a Boeing 757 from an alti- 
tude of about 1 2,000 ft, where it exits the TRACON, to 34,000 
ft, in a randomly selected wind field. The solid line represents 
the best-fit parabola, and the dashed lines represent an 
15 example error tolerance of ±2000 ft relative to the parabola. 
The constant-CAS segment is followed by a short constant - 
Mach segment (not shown), which would require its own 
curve fit. In most cases, the aircraft should be able to fly the 
specified CAS of 296 knots without altitude feedback or 
20 throttle modulation, and stay within the specified altitude 
range. Only if the TS is substantially in error would the 
aircraft need to use feedback of altitude, and perhaps also 
throttle modulation, to stay within tolerance. Such error could 
be due to errors in wind, thrust, and/or weight. Altitude feed- 
25 back, and perhaps throttle modulation, could be activated 
when the altitude deviation from the reference reaches a 
threshold value of, say, 1 000 ft. Alternatively, an FMS could 
be programmed to use a more refined criterion. 

The curve fit error bounds of the parabola in FIG. 4 are 
30 -1 89 ft and +289 ft, for atotal range of 478 ft; this curve could 
be offset to make the error bounds symmetric, if desired but 
that is not done here). With an ALT error tolerance of ±2000 
ft. that fit allows a worst -case altitude deviation, relative to the 
TS output, of 189-2000=-181 1 to2000-289=+171 ft, which 
35 is probably sufficient. However, if the error tolerance were 
tighter, say ±1000 ft, that quadratic fit would only leave a 
worst-case altitude deviation of -811 ft to +711 ft, which 
might not be considered sufficient. In that case, the segment 
could be divided into two or more segments, or a cubic or 
40 quartic polynomial could be used for a better fit. For this 
example, a cubic polynomial gives fit error bounds of -178 to 
+94 ft (27 ft range), and a quartic gives -102 to +68 ft (170 ft 
range). Polynomials of fifth order or higher may have numeri- 
cal problems and should be avoided, but polynomials of 
45 fourth order or less should not suffer from significant numeri- 
cal round-off errors if a consistent numerical precision of 64 
bits is used in both ground-based and airborne computers. 
When engine problems prevent climbing at a normal rate, the 
aircraft should notify the ground immediately of the expected 
50 non-conformance, and the ground will then reroute any 
affected aircraft. When routing under climbing aircraft, pre- 
cautions could be taken to minimize the chance of a conflict 
due to engine non-performance. For example, trajectories 
could be determined for all levels of engine power from full 
55 power down to one engine out and those trajectories could be 
avoided. The lower altitude tolerance, which will be dis- 
cussed later, could be used for this purpose. 

FIG. 5 shows the AT location associated with the climb of 
FIG. 4. The best-fit parabola fits with error bounds of -0.05 to 
60 +0.07 nmi, which is close enough for all practice purposes. 
The example error tolerances represented by the dashed lines 
start out at ±2 nmi and grow linearly with time at a rate of 0.5 
nmi/min to ±6 nmi at 8 min from the start of the climb. 

FIG. 6 shows the ALT profile synthesized by the TS for a 
65 constant-CAS, idle-thrust descent segment of a Boeing 727 
from an altitude of about 30,000 ft to 1 1 ,000 ft, where it enters 
the TRACON, in a randomly selected wind field. Again, the 
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solid line represents the best- fit parabola, and the dashed lines 
represent a hypothetical error tolerance of ±1500 ft. The 
constant CAS-segment is preceded by a short constant-Mach 
segment (not shown), which would require its own curve fit. 
Again, the aircraft should normally be able to fly the constant 
CAS of 280 knots without altitude feedback or throttle modu- 
lation and stay within the specified altitude range. As before, 
altitude feedback, and perhaps throttle modulation, could be 
activated when the altitude deviation reaches some threshold 
value. With error bounds of -53 ft to +1 1 6 ft, the curve fit for 
this descent is much more accurate than for the climb of FIG. 
4. Descents tend to be more nearly linear than long climbs, 
and are usually well modeled with a parabola. In general, an 
arrival descent would be followed by a short level cruise 
segment into the meter fix, which would allow the aircraft to 
cross the meter fix at a precise level altitude. 

FIG. 7 shows the AT location associated with the descent of 
FIG. 6. The best-fit parabola fits with error bounds of -0.06 to 
+0.06 nmi, which is close enough for all practical purposes. 
The example error tolerances represented by the dashed lines 
start out a 6 nmi and are reduced in magnitude linearly with 
time at a rate of 0.5 nmi/min to slightly less than ±2 nmi at 9 
min from the start of descent. Climbing aircraft usually only 
need to control airspeed using the elevator (with throttle 
fixed). However, descending aircraft neec to control both 
altitude and precise time of arrival at the meter fix, so that such 
aircraft need tighter control using not only the elevator but 
possibly also the throttle or speed brake. 

A wind vector is added to the airspeed “vector” (based on 
airspeed and heading relative to the wind) to determine the 
varying groundspeed of the aircraft as it progresses through a 
turn. The law of cosines is then used to compute the resulting 
groundspeed, which can be numerically integrated to pre- 
cisely determine the AT location as a function on time. In any 
nonzero wind field the groundspeed will vary through the 
turn, and the bank angle must be varied accordingly to main- 
tain a coordinated turn (in which gravity and centrifugal 
acceleration add vectorially to a force normal to the floor of 
the aircraft). The necessary bank angle (3 is given in terms of 
the varying groundspeed, v, by 

P=tan -^/(Rjg)), (9) 

where R r is the (constant) turn radius, and g is gravitational 
acceleration. The radius of the turn can be selected so that the 
maximum bank angle through the turn does not exceed a 
specified magnitude. 

The turn angle between adjacent straight (great circle) 
segments can be determined using the “course” function 
given in the Appendix for determining the initial course angle 
of a great circle segment (the spherical Earth equations are 
plenty accurate for this purpose The turn angle between great 
circles (A,B) and (B,C) is “course(B,C)-course(B,A)±180 
deg.” Each turn segment should be tangent to the two great 
circle segments that it connects The turn corner is at the 
intersection of the two great circle segments, and a turn of 
angle A must begin at a distance d=r tan(A%) before the 
corner, as shown in the FIG. 8 turn geometry. 

A problem with transitioning from a non-level segment to 
a level segment is that the altitude tolerance will virtually 
always be discontinuous. A typical altitude tolerance for level 
flight might be ±200 ft. but for climb or descent the tolerance 
could be ten times larger. Going from a level segment to a 
non-level segment is not a problem because the tolerance 
increases, but going from a non-level segment to a level 
segment is a problem because the altitude tolerances decrease 
sharply and discontinuously, as illustrated in FIG. 9. The 
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linearly decreasing “transition tolerance” shown for the lower 
altitude tolerance is one possible approach for reducing the 
tolerance less abruptly. The AT distance over which the tol- 
erance decreases, which could be something like 5 to 10 nmi, 
5 will be specified as part of the segment characteristics for the 
level segment. 

FIG. 9 also illustrates another problem with transitioning 
from a non-level segment to a level segment. The upper alti- 
tude tolerance during the climb segment allows the aircraft to 
to go significantly above its intended cruising altitude, exposing 
it to potential conflicts with traffic at the next higher flight 
level. To prevent such exposure, an “altitude limit” can be 
specified in the “max” attribute of the “ALT” sub-element of 
“segment.” This maximum altitude limit overrides the upper 
15 altitude tolerance, as shown in the figure. For a descent seg- 
ment, the “max” attribute would be replaced by “min.” This 
maximum altitude limits any “overshooting” of target altitude 
and will normally be the upper altitude tolerance for the level 
segment. The same geometry turned upside down applies to 
20 descent. 

The three axes in which the error tolerances are specified 
are AT, CT and ALT (vertical). In terms of energy, the least 
expensive axis in which to maintain conformance is the cross - 
track axis. If the aircraft can fly the assigned ground track to 
25 within the CT tolerance, the marginal energy required to do so 
is usually small. In the case of unexpected crosswinds, the 
aircraft only has to “crab” to stay on track, and the main 
controls required are the ailerons and the rudder. 

In a cruise mode, vertical conformance is also inexpensive, 
30 but AT conformance can be far more expensive, depending on 
how accurately the AT winds can be modeled and predicted. 
Aircraft normally cruise at constant airspeed (CAS or Mach), 
and the groundspeed corresponding to the most efficient air- 
speed obviously varies with the AT wind speed. If the wind 
35 field prediction is accurate, then an efficient groundspeed can 
be determined and assigned. If the predicted winds are sub- 
stantially in error, however, the airspeed corresponding to the 
assigned groundspeed could be inefficient or even realistic. 
The two relevant concerns here are the wind prediction accu- 
40 racy and the speed range of the aircraft. 

R. E. Cole, Green, Jardin and Schwartz (“Wing Prediction 
Accuracy for Air Traffic Management Decision Support 
Tools, Third USA/Europe Air Traffic Management R&D 
Seminar, Napoli, Italy, 13-16 Jun. 2000) cite wind prediction 
45 accuracy results for the Rapid Update Cycle (RUC-1) aug- 
mented with aircraft wind reports. The wind error vector 
magnitude was 7.85 mis (15.3 knots) or less 90 percent of the 
time, and error vector magnitude exceeded 10 mis (19.4 
knots) only 4 percent of the time. The errors tend to be 
50 somewhat worse during the winter months because of higher 
wind speeds in general, but they are somewhat better than the 
quoted figures the rest of the year. The errors also tend to be 
larger at higher altitudes where the wind speeds are higher, 
but the figures quoted above are from actual aircraft at their 
55 operating altitudes. These performance figures can perhaps 
be expected to improve over the next twenty years. Note also 
that the AT component of the wind error, which is the signifi- 
cant quantity here, is less than the vector magnitude. Aver- 
aged over all heading directions, the mean headwind error is 
60 2/,rr('; : : , 0.64) times the magnitude of the error vector (that 
probably overstates the average effective reduction, however, 
because neither heading directions nor wind error vectors are 
uniformly distributed). 

The other relevant factor here is aircraft speed range. FIG. 
65 1 0 is a graph of the speed envelope for an MD-80 aircraft at a 
gross weight of 135,000 lb on a standard day. This Figure is 
based on data that came from a McDonnell Douglas perfor- 
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mance handbook. These data, as well as the aircraft weight, 
would be known by the ground systems for all equipped 
aircraft. Any uncertainties in the weight would need to be 
accounted for to guarantee that the speed range is not over- 
estimated. 

FIG. 10 indicates that the speed range decreases sharply as 
the pressure altitude ceiling of approximately FL330 (for this 
weight) is approached. At FL310 the speed range is approxi- 
mately 75 knots, and at FL290 it is slightly below 100 knots. 
Assume that the aircraft is flying at FL310 at the recom- 
mended cruise speed of Mach 0.76, which is equivalent to 446 
knots at that altitude. The speed envelope then goes from 
approximately 394 to 469 knots, as shown in FIG. 10 , so the 
speed can be increased by a maximum of 23 knots or 
decreased by a maximum of 52 knots from the recommended 
speed. Thus, if the wind prediction error is within the range of 
-23 to +52 knots, the aircraft can maintain the assigned 
groundspeed exactly. 

Although the aircraft can fly from 394 to 469 knots at 
FL310, it cannot fly efficiently over that entire range, of 
course. Suppose the efficiency is deemed “acceptable” from 
430 to 454 knots at that altitude and weight. Then, as long as 
the wind prediction error is within -8 to + 1 6 knots, the aircraft 
can maintain the assigned groundspeed exactly and still fly 
with acceptable efficiency. Outside that range, an incentive 
exists to fly at 500 an efficient airspeed until a point is reached 
at which the aircraft can no longer conform to the AT error 
tolerances. Hence, additional rules may need to be estab- 
lished for how tightly an aircraft should track the assigned 
groundspeed, but such rules will not be discussed in this 
disclosure. 

The aircraft isn’t required to fly the assigned groundspeed 
exactly; it is only required to stay within the AT location 
bounds, which can vary linearly with time. Suppose the AT 
tolerance rates are zero, the AT location tolerances are ±2 nmi, 
and the aircraft is centered within the AT location bounds. If 
the AT wind error is within ±10 knots, the aircraft will be able 
to maintain its recommended airspeed for at least 12 min, 
worst case, before it can possibly fall out of conformance. 

Alternatively, if the along-track tolerance rates are ±12 
knots, then the AT tolerances will increase at that rate, which 
is 0.2 nmi per minute, in each direction. Thus, the aircraft will 
be able to fly at the recommended airspeed and maintain a 
constant or increasing margin from the AT bounds even if the 
AT wind predictions are in error by up to ±12 knots. 

Commercial transport airplanes normally climb with 
throttle fixed somewhere in the range of 85 to 95 percent of 
full throttle, with feedback to the elevator to maintain con- 
stant CAS (at lower altitudes) or constant Mach (at higher 
altitudes). If the errors in the prediction of the AT wind speed 
are reasonably small, the aircraft should still be able to fly in 
that mode. However, when ALT or AT location approach their 
bounds, the throttle may also need to be adjusted to maintain 
conformance. The feedback to the throttle could be pro- 
grammed to start automatically when the deviation reaches 
some threshold magnitude. The error tolerances should be set 
to accommodate the entire range of potential wind errors to 
some high level of certainty, say 99.9 percent. This could 
make the error tolerances fairly large, but at least they will 
clearly bound the area (as a function of time) that needs to be 
avoided by other aircraft. In the current system, the lack of 
explicit bounds forces controllers to effectively block out 
excessively large amounts of airspace for climbing and 
descending aircraft, which reduces airspace capacity. 

In the event of substantial errors in the prediction of AT 
winds, aircraft could be forced to fly at grossly inefficient 
speeds to maintain AT location conformance. Worse yet, they 
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could reach a state in which they are aerodynamically inca- 
pable of flying at the speed necessary to maintain conform- 
ance. To avoid either of those two undesirable conditions, 
particularly the latter, the AT assignments will be updated 
5 periodically. The updates could apply to position, speed, and 
error tolerances. With proper updates, the worst that should 
happen is that some traffic may be forced to fly inefficiently 
for short periods of time to avoid a conflict, but they would 
obviously never be required to fly at speeds of which they are 
to incapable of flying. 

A complete 4D trajectory can (optionally) be filed by each 
participating aircraft or its AGC prior to takeoff. The trajec- 
tory reference time can be defined relative to a time at which 
the aircraft is expected to cross some predefined marker, such 
1 5 as the end of the takeoff runway. Because takeoff time usually 
cannot be predicted exactly, the first AT update will occur 
immediately after takeoff. When the aircraft crosses the ref- 
erence marker, its reference time will be adjusted accord- 
ingly. Because all other times are relative to the reference 
20 time, no other times need to be changed. By adjusting the 
trajectory reference time, the entire trajectory can be effec- 
tively shifted in time. 

After takeoff, the FMS will guide the aircraft along its 
assigned climb trajectory. As explained earlier, conventional 
25 feedback of speed error to the elevator will be used to main- 
tain constant CAS or Mach, and feedback to the throttle will 
be used only if the vertical or AT deviation reaches some 
threshold value, which shouldn’t happen often if the wind 
predictions are reasonably accurate and the error tolerances 
30 are reasonable. If the aircraft does drift away from its refer- 
ence trajectory and approach its AT error bounds, however, 
the ground can update the AT assignment by changing the 
assigned position, speed, and/or error tolerances. Such 
updates would be done only if they do not cause a conflict 
35 within the conflict time horizon of, say, 1 5 min. 

The most common type of AT assignment update will be to 
change the assigned position to the current position and to 
simultaneously reset the position error tolerances to their 
initial values. The assigned groundspeed (the rate of change 
40 of the assigned AT location) could also be changed if the wind 
model is determined to be significantly in error. This kind of 
update could be done periodically at a rate of, say, once per 
two minutes, except that it would not be done if it produces a 
potential conflict within the conflict time horizon. A potential 
45 conflict is defined as having the bounding spaces of two 
aircraft come closer together than the required minimum 
separation. In other words, conformance to their assigned 
trajectories by any two aircraft must guarantee the minimum 
required separation. 

50 Because all the times given in the trajectory specification 
are relative to the trajectory reference time, the entire trajec- 
tory can be shifted in time by changing the reference time. 
That is equivalent to changing the assigned AT location. As a 
flight progresses, trajectory segments that are completely in 
55 the past can be discarded. Because AT updates will be a 
common operation, an abridged format is appropriate, where 
unchanged data is not repeated. 

Trajectory Visualization. 

60 As for the unique onboard requirements, the main one is 
that the aircraft be capable of conforming to its assigned 
trajectory. This will most likely be an automated function of 
an advanced FMS. Another onboard requirement is to display 
the reference trajectory for the pilot, which could be done 
65 with an advanced CDTI (Cockpit Display of Traffic Informa- 
tion), perhaps as a head-up display. A CDTI could also con- 
ceivably be used for real-time guidance, instead of an FMS, to 
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allow aircraft not equipped with an FMS to fly in high-density 
airspace for short periods of time for arrival at busy airports. 
This possibility will be discussed briefly in the next section. 

Trajectory Display 

In twenty years, virtually all commercial transport aircraft 
are likely to be equipped with an advanced FMS capable of 
automatically flying a specified trajectory. Pilots will still 
need to be able to visualize their assigned trajectory and 
monitor conformance, however, New. CNS technologies 
should be able to provide that capability at a reasonable cost. 

A technology called Cockpit Display of Traffic Information 
(CDTI) uses GPS/WAAS and ADS-B to provide an onboard 
visual representation of the local traffic environment. An 
advanced CDTI could also show assigned trajectories for 
improved situational awareness. 

Such an advanced CDTI could also conceivably be used in 
place of an FMS for pilot guidance. Whereas an FMS must be 
developed for a particular aircraft model, generic CDTI units 
can be produced for any aircraft with a compatible slot in the 2Q 
cockpit, and are potentially cheaper to test and certify, hence 
a CDTI unit could potentially be cheaper than an FMS. This 
fact could allow lower-cost aircraft, perhaps even general 
aviation, to participate for short periods of time as equipped 
aircraft in high-density airspace, which could be critical in 25 
crowded terminal and transition airspace. The trajectory error 
tolerances for such aircraft would probably need to be larger 
than for aircraft equipped with an FMS, but they would still 
fly much more precise and predictable trajectories than air- 
craft not equipped with CDTI guidance. 30 

FIG. 11 shows how an advanced CDTI might be used for 
situational awareness. It shows the local traffic as any CDTI 
would, but it also shows the bounding space and the future 
assigned horizontal path of the user and the local traffic, so 
surprise maneuvers are minimized. The start of descent point 35 
is also indicated, and a vertical guidance bar appears on the 
right edge of the display so that all three axes are represented 
in one view. The small box in the upper-left corner shows 
current altitude, course, and speed, and a compass pointer 
appears in the lower-left comer. 40 

FIG. 12 shows how the same CDTI could be zoomed in for 
finer guidance, again perhaps as a head-up display. The hori- 
zontal position of the user within the bounding space is 
clearly shown, as is the velocity relative to the reference 
velocity. Again, a vertical guidance bar appears on the right 45 
edge of the display, but more detailed vertical guidance could 
also be provided during climb and descent, as shown in FIG. 

13. The vertical guidance display includes a CT guidance bar 
so that some guidance is always provided in all three axes. 

Each AOC or pilot will be allowed to submit requests at any 50 
time for specific trajectories or trajectory revisions. 
Requested trajectories and revisions will be checked on the 
ground for conflicts within the conflict time horizon, and if 
they are free of conflicts and consistent with the longer range 
traffic flow plan, they will be approved. Otherwise, the 55 
requested trajectory will be minimally revised on the ground 
to resolve the conflict or other problem, then it will be 
uplinked as the assigned trajectory. Message digests such as 
MD5 or SHA-1 can be used to guarantee datalink integrity. 

Trajectories are broken down into a series of trajectory 60 
segments of various types. With three vertical types (climb, 
level, and descent), three heading types (straight, right turn, 
and left turn), and five speed types, the total number of pos- 
sible combinations is 45. The horizontal path consists of a 
series of geodetic waypoints connected by great circles, and 65 
the great-circle segments are connected by turns of specified 
radius. Along -track location is specified as a low-order poly- 
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nomial function of time, and vertical profiles for climb and 
descent are specified as low-order polynomial functions of 
actual (not reference) AT location. Flight technical error tol- 
erances in the AT, CT and ALT axes determine the bounding 
space. Periodic updates in the along-track axis adjust for 
errors in the predicted AT winds. 

Use of this regimen of assigned 4D trajectories can elimi- 
nate the need for separation monitoring (of equipped aircraft 
by human air traffic controllers. It can also guarantee that the 
equipped traffic will be able to fly free of conflicts for at least 
several minutes even if all ground systems and the entire 
communication infrastructure fail. This failsafe guarantee, 
along with the elimination of the human factor from the 
primary separation feedback loop, has the potential to greatly 
increase airspace capacity. That increase in capacity will be 
necessary within the next two decades if the U.S. air traffic 
system is to keep pace with the growing demand for air travel . 

The invention claimed is: 

1 . A system for analyzing and presenting information on an 
aircraft flight path, the system comprising a computer that is 
programmed: 

to provide a description of aircraft location along an actual 
flight path in terms of Cartesian coordinates (x, y, z), 
referenced to a selected coordinate system, and in terms 
of an estimated distance 6 moved by the aircraft along a 
reference flight path during an elapsed time interval 
corresponding to the present time; 
to estimate coordinates (x , y p , z ) for present location for 
the aircraft; 

to estimate location of a perpendicular foot, having loca- 
tion coordinates (x^ 00? , y foot ^ z /bo J of the aircraft on a 
curve C representing the reference flight path; 
to interpret the perpendicular foot location as a nearest 
reference flight path location corresponding to the 
present location of the aircraft; 
to estimate an arc length along the curve C from an initial 
reference location to the nearest flight path location; 
to interpret the arc length as the estimated distance 5; 
to estimate, at each of a sequence of values of the estimated 
distance 6, an error value, including at least one of an 
along-track (“AT”) error, a cross-track (“CT”) error and 
an altitude (“ALT”) error for the aircraft relative to the 
reference flight path; 

to provide, for at least one of the sequence of values of 5, a 
permitted error range for at least one of the AT error, the 
CT error and the ALT error, wherein the AT error, the CT 
error and the ALT error define a bounding space that 
maintains a minimum required separation distance; and 
to present at least one of a visually perceptible display, an 
audibly perceptible display and an alphanumeric display 
of at least one of the AT error, the CT error and the ALT 
error as an error value for at least two of the sequence of 
values of 6. 

2. The system of claim 1, wherein said computer is further 
programmed: 

to advise an operator of said aircraft of a magnitude of said 
error value; and 

to recommend an action to be taken by said aircraft to 
correct or reduce the magnitude of said error value, 
when said error value lies outside a corresponding per- 
mitted error range. 

3. The system of claim 1, wherein each of M locations 
(M^2), numbered m=l, . . . , M, on said curve C qualifies as 
said perpendicular foot for said aircraft location on said curve 
C with a corresponding estimated distance §(foot;m), and 
said computer is ftirther programmed to estimate a designated 
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perpendicular foot, for purposes of estimating said arc length 
along said curve C, by a process comprising: 

(1) estimating said distance 8(foot;m) (m=l, . . . , M) for the 
perpendicular foot on said curve C for the location num- 
ber m; and 

(2) the designated perpendicular foot is determined for a 
location number m ml (l^ml^M) for which 6(foot; 
ml) is no greater than any of the distances 6(foot;m) 
(m=l, . . . , M). 

4. The system of claim 1, wherein said computer is further 
programmed: 

to estimate a second distance S(PLAET) that said aircraft 
would have flown along said reference flight path in said 
elapsed time interval; and 

to estimate said AT error as a distance 8(AT)=l8-8- 
(PLAET)I. 

5. The system of claim 1, wherein said computer is further 
programmed: 

to estimate a second distance 8(PLAET), having location 
coordinates i^ PLAEP Y PE aet> z plaet)> that said aircraft 
would have flown along said reference flight path in said 
elapsed time interval; and 

to estimate said CT error as a distance ^(CT^dllKx^y^)- 
(^PLAETtf plaet)% where dl || || is a selected distance met- 
ric for two locations in a plane having location coordi- 
nates (x p ,y p ) and plaet)- 

6 . The system of claim 1 , wherein said computer is further 
programmed: 
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to estimate a second distance S(PLAET), calculated from 
location coordinates i^ PLAEP y PLAEP z plaet)> that said 
aircraft would have flown along said reference flight 
path in said elapsed time interval; and 
5 to estimate said ALT error as a distance 8(ALT)=l(z ), 
( z plaet)\- 

7. The system of claim 1, wherein said computer is further 
programmed: 

to estimate a second distance 8(PLAET), having location 
10 coordinates (x PLAET , y PLAET , z PE aei)> that said aircraft 
would have flown along said reference flight path in said 
elapsed time interval; 

to estimate said AT error as a distance 8(AT)=l8-8 
(PLAET); 

15 to estimate said CT error as a distance 8(CT) dl||(x ,y )- 
( x PLAEpy plaet) 1 1’ where dl|| || is a selected distance met- 
ric for two locations in a plane having location coordi- 
nates (x p ,y p ) and (x PLAET ,y PLAET ); and 

to estimate a horizontal error as a distance 

20 8(horiz)={8(AT) 2 +8(CT) 2 } 1/2 . 

8. The system of claim 7, wherein said computer is further 
programmed: 

to estimate said ALT error as a distance 6(ALT)=l(z /5 ), 
25 ( z plaet) I j an d 

8(total)={8(AT) 2 +8(CT) 2 +8(ALT) 2 } 1/2 . 



